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C _^ Abstract. The main objective of this article is to study the nature of the 

£N| Andrews critical point in the gas-liquid transition in a physical-vapor trans- 

I port (PVT) system. A dynamical model, consistent with the van der Waals 

^ equation near the Andrews critical point, is derived. With this model, we de- 

K - 5 duce two physical parameters, which interact exactly at the Andrews critical 

point, and which dictate the dynamic transition behavior near the Andrews 
critical point. In particular, it is shown that 1) the Andrews critical point 
is a switching point where the phase transition changes from the first order 
to the third order, 2) the gas-liquid co-existence curve can be extended be- 
yond the Andrews critical point, and 3) the liquid-gas phase transition going 
beyond Andrews point is of the third order. This clearly explains why it is 
ry^ hard to observe the gas-liquid phase transition beyond the Andrews critical 

. t point. Furthermore, the analysis leads naturally the introduction of a general 

C^ asymmetry principle of fluctuations and the preferred transition mechanism 

C for a thermodynamic system. 
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1. Introduction 

Phase transition is one of the central problems in nonlinear sciences. Many 
fvi systems have different phases, and the most commonly encountered phases are 

£SJ gas, liquid and solid phases. A natural system which possesses these three phases 

is the physical-vapor transport (PVT) system. As we know, a PVT system is a 

system composed of one type of molecules, and the interaction between molecules is 

(^. governed by the van der Waals law. The molecules generally have a repulsive core 

and a short-range attraction region outside the core. Such systems have a number 
of phases: gas, liquid and solid, and a solid can appear in a few phases. The most 
typical example of a PVT system is water. 

A PT-phase diagram of a typical PVT system is schematically illustrated by 
Figure |1.1| where point A is the triple point at which the gas, liquid, and solid 
;Jj phases can coexist. Point C is the Andrews critical point at which the gas- liquid 

coexistence curve terminates [3l 0] . Classical view on the termination of the gas- 
liquid coexistence curve at the critical point amounts to saying that the system 
can go continuously from a gaseous state to a liquid state without ever meeting an 
observable phase transition, if we choose the right path. 
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Figure 1.1. Coexistence curves of a typical PVT system: A is 
the triple point, C is the critical point, and the dashed cure is a 
melting line with negative slope. 



It is, however, still an open question why the Andrews critical point exists and 
what is the order of transition going beyond this critical point. In pQ, a mathe- 
matical theory is derived to address this problem. In this article, we explore the 
physical implications of the mathematical theory derived in pQ, give a theory on 
the nature of the Andrews critical point, and introduce the asymmetry principle of 
fluctuations and the preferred transition mechanism. 

First, the modeling is based on 1) the Landau mean field theory, 2) a unified 
dynamic approach for equilibrium phase transitions, 3) the classical phase diagram 
in Figure |1.1| and c) the van der Waals equation. It is worth mentioning two im- 
portant aspects of the model we derived. One is that the new model can be used to 
study liquid-solid and gas-solid transitions as well, by choosing different parameters. 
Second is the consistency of the model with the van der Waals equation. Namely, 
near the Andrews critical point C, the steady state equation of the homogeneous 



model (2.12) is exactly the van der Waals equation. This consistency gives a good 



validation of the mean field model. In addition, the dynamic approach leading to 



the model provides much richer information. For example, the model (2.7) can be 
used to study the heterogeneity of the system. 

Second, with the dynamic model at our disposal, we introduce two new physical 
parameters A = X(T,p) and a 2 = a,2(T,p), where the temperature T and the 
pressure p are control parameters. These two physical parameters determine the 
phase transition behavior near the Andrews point and provide the key ingredient 
to characterize the nature of the Andrews point. 

It is remarkable that these two parameters reproduce the location of the Andrews 
critical point C, which is the same as derived by van der Waals in his classical work, 
although the method we use is the dynamic approach based on the Landau mean 
field theory, different from the one used by van der Waals. Coincident ally, these 
two parameters are the second and third-order derivatives of the Gibbs energy at 



the equilibrium state po] see (3.3) 



Third, the two parameters A and a 2 correspond to two-curves in the pT-phase 
plane, and interact exactly at the Andrews critical point C '. Then with the dynamic 
transition theory developed recently by the authors, we deduce a theory on the 
Andrews critical point C: 1) the critical point is a switching point where the phase 
transition changes from the first order with latent heat to the third order, and 2) the 
gas-liquid phase transition beyond the Andrews critical point is of the third order. 
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This explains why it is hard to observe the phase liquid-gas transition beyond the 
Andrews point, and clearly 

Fourth, physical intuition and the theory lead us to introduce asymmetry prin- 
ciple of fluctuations, and the preferred transition mechanism. 

2. A Dynamic Model for Gas-Liquid Transition 

The classical and the simplest equation of state which can exhibit many of the 
essential features of the gas-liquid phase transition is the van der Waals equation: 

(2.1) (p+ ^){v-b)=RT, 

where v is the molar volume, p is the pressure, T is the temperature, R is the 
universal gas constant, b is the revised constant of inherent volume, and a is the 
revised constant of attractive force between molecules. If we adopt the molar density 



p = 1/v to replace v in (2.1), then the van der Waals equation becomes 
(2.2) - (bp + RT)p + ap 2 - abp 3 + p = 0. 

Now, we shall apply thermodynamic potentials to investigate the gas-liquid phase 
transitions in PVT systems, and we shall see later that the van der Waals equation 
can be derived as a Euler-Langrange equation for the minimizers of the Gibbs free 
energy for PVT systems at gaseous states. 

Consider an isothermal-isopiestic process. The thermodynamic potential is taken 
to be the Gibbs free energy. In this case, the order parameters are the molar density 
p and the entropy density S, and the control parameters are the pressure p and 
temperature T. The general form of the Gibbs free energy for PVT systems is 
given as 



dx. 



(2.3) G(p,S,T,p) = J [^\Vp\ 2 + ^\VS\ 2 +g(p 1 S 1 T 1 p)-ST-a(p,T,p)p 

where g and a are differentiable with respect to p and S, Q C K 3 is the container, 
and ap is the mechanical coupling term in the Gibbs free energy, which can be 
expressed by 

(2.4) a(p, T,p)p = pp- ^bp 2 p, 

where b = b(T,p) depends continuously on T and p. In fact, this mechanical 



coupling term should be p. In view of the van der Waals equation (2.2) and the 



mathematical analysis based on the new dynamical transition theory, phenomeno- 



logically we need to adjust the term by adding a coefficient a, leading to (2.4) as 
the first two terms in the Taylor expansion. Although the van der Waals equation 
works for gaseous sates only, by choosing the dependence of the coefficient b on the 
temperature and the pressure, the energy applies to the liquid and solid states as 
well. This is a very subtle term from the physical point of view to derive a feasible 
free energy. 

Based on both the physical and mathematical considerations, we take the Taylor 
expansion of g{p, S, T,p) on p and S as follows 

(2.5) g = l -a lP 2 + ^S 2 + (3 2 Sp 2 - ^a 2 p 3 + *a 3 p 4 , 
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where on (1 < i < 3), /?i and /?2 depend continuously on T and p, and 

(2.6) f3 1 =p 1 (T,p)>0, a l = a l (T,p)>0 » = 2,3. 

In a PVT system, the order parameter is u = (/?, 5) , 

p = p\ — pa, S = Si — So, 

where pi and Si (i = 0,1) represent the density and entropy, Po,Sq are reference 
points near the coexistence curve of gas and liquid states. Hence the conjugate 
variables of p and S are the pressure p and the temperature T. Thus, by the le 



Chatlier principle, we derive from (2.3|-(2.5| the following dynamic model for a 
PVT system: 



(2.7) 



dp 

dt 

dS_ 



piAp - (ai + bp)p + a 2 p 2 - a 3 p 3 - 2(3 2 pS + p, 
-p 2 AS-(3 1 S-p2P 2 +T. 



A physically meaningful boundary condition for the system is the Neumann 
boundary condition: 



(2.8) 



d P n dS n an 

= 0, ^— =0 on oil. 



dn dn 

An important special case for PVT systems is that the pressure and temperature 
functions are homogeneous in 51. Thus we can assume that p and S are independent 



of x £ n, and the free energy (2.3) with (2.4) and (2.5) can be expressed as 

bp 2 p 



(2.9) G( P ,S,T,p) 



"12 



ft rf 



S' + p 2 s P z 



Y p3 + T p4 



pp - ST. 



From (2.9) we get the dynamical equations as 
(2.10) 



dp 
~di 
dS_ 
~dt 



(ai + bp)p + a 2 p 2 - a 3 p 3 



2f3 2 S P + p, 



-PiS-(3 2P 2 +T. 



Because /3i > for all T and p, we can replace the second equation of (2.10) by 

5 = /?r l (T-/? 2 p 2 ). 



(2.11) 



Then, (|2. 10|) are equivalent to the following equation 

(ai +bp + 2^ 1 p 2 T)p + a 2P 2 - (a 3 - 2^ 2 /3^)p 3 +p. 



(2.12) 



dp 
~df 



It is clear that if a\ — 0,2/3 1 1 (3 2 = R,a 2 — a, (03 — 2/3f/3i 1 ) = ab, then the 
steady state equation of (2.12) is referred to the van der Waals equation. 



We remark that (2.12) can be considered as the dynamic version of the van der 



Waals equation, although we used the Landau mean field theory together with 
the le Chatlier principle. The approach provides a much richer information. For 



example, the model (2.7) can be used to study the heterogeneity of the system. In 



addition, the model here can be used to study liquid-solid and gas-solid transitions 
as well, by choosing different parameters. 
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3. Two New Physical Parameters and the Andrews Critical Point 



In this section we use (2.12| to derive two new physical parameters, which dic- 



tates the dynamic transition behavior near the Andrews critical points. 



Let po be a steady state solution of (2.12| near the Andrews point C = (T c ,p c ). 
We take the transformation 

P = Po + p'- 



Then equation (2.12) becomes (drop the prime) 



(3.1) — = Xp + a 2 p 2 - abp 3 . 

where 

A = 2apo — 3abpl — a\ — bp — RT, 
a 2 = a(l - 36p ). 

where a% is close to zero. Here we emphasize that po & n d (A, a 2 ) are all functions 
of the control parameters (T, p) . 

These are two important physical parameters, which are used to fully characterize 
the dynamic behavior of gas-liquid transition near the Andrews point. In fact, from 
the derivation of the model, we obtain immediately the following physical meaning 
of these two parameters: 

( 3 - 3 ) X ^ = dp^Uo' fl2(T ' P) = 2^l.=Po' 

where po = Pa(T,p) is the equilibrium state. 

In the PT-plane, near the Andrews point C = (T c ,p c ), the critical parameter 
equation 

X = X(T,p) = in \T-T e \<8, \p- Pc \<5 
for some 5 > 0, defines a continuous function T = <f>(p), such that 

<0 \lT>cj){p), 

(3.4) \[ =0 if T = ^(p), 

>0 iiT<(j)(p). 

Equivalently, this is called the principle of exchange of stabilities, which, as we 
have shown in iT], is the necessary and sufficient condition for the gas-liquid phase 
transition. 

One important component of our theory is that the Andrews critical point is 
determined by the system of equations 

A = 0, 

(3.5) «2 = 0, 

- (bp + RT)p + apl - abpl + p = 0. 

Here the first equation is the critical parameter equation, the second equation, as we 
shall see below, determines the switching point where the phase transition switches 
types, and the last equation is the van der Waals equation, which is also the steady 
state equation of the dynamic model. 
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Figure 4.1. The point C = 
(T c ,p c ) is the Andrews critical 
point. 




Then by a direct computation, it is easy to see that the critical point C is given 



by 

(3.6) 



(p c ,T c ,p c 



la 8a \ 
36'2762'276i?J 



This is in agreement with the classical work by van der Waals. Here we obtain the 
Andrews point using a dynamic approach. 

4. Theory of the Andrews Critical Point 
We now explain the gas-liquid transition near the Andrews critical point C. 



First, we have shown in ( |3.6[ ) that at the equilibrium point po, the two curves 
given by \(p,T) = and ci2{p,T) = interact exactly at the critical point C as 
shown in Figure |4.1| and the curve segment AB of A = is divided into two parts 
AC and CB by the point C such that 

a 2 (T,p)>0 for (T,p) £ AC, 

a 2 (T,p)<0 for (T,p)eCB. 

Here the curve AC is the classical gas-liquid co-existence curve. 

Second, on the curve AC, excluding the critical point C, a 2 (T,p) > 0. The 
phase transition of the system is a mixed type if we take a pass crossing the AC; 



see Figure 4.2 In Figure 4.2 p is the deviation from the basic gaseous state po- 
We now consider different states given in Figure |4.2|[a) : 




Figure 4.2. Type-Ill (mixed) transition for a 2 > 0: (a) Phase 
diagram for fixed pressure, and (b) phase diagram for fixed tem- 
perature. 



(1) For T > Ti, the gaseous state po, corresponding to p = in the figure, is 
stable. This is the only stable physical state in this temperature range, and 
the system is in the gaseous state. 
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(2) For To < T < Ti, there are two metastable states given by gaseous phase 
Po and the liquid phase p + p + . 

(3) For T < T , there are three states: the unstable basic gaseous state po, 
and the two metastable states: po + p~ and po + P + ■ One important 
component of our theory is that the only physical phase here is the liquid 
phase represented by metastable state: po + p + . Although, mathematically 
speaking, the gaseous state po + p~ is also metastable, it does not appear in 
nature. The only possible explanation for this exclusion is the asymmetry 
principle of fluctuations, to be further explored in the next section. 

(4) Hence we have shown that as we lower the temperature, the system un- 
dergoes a first order transition from a gaseous state to a liquid state with 
an abrupt change in density. In fact, there is an energy gap between the 
gaseous and liquid states: 

AE = G( Po + p+) G( P o) ~ -\{ P + f !|(p + ) 3 < for T < Ti; 

see pp. This energy gap \AE\ stands for a latent heat, and AE < shows 
that the transition from a gaseous state to a liquid state is an isothermal 
exothermal process, and from a liquid state to gaseous state is an isothermal 
cndothermal process. 

Third, at the critical point C, we have 02 = 0. Then the dynamic transition 



is as shown in Figure 4.3 see [T] for the detailed mathematical analysis leading to 
this phase diagram: 

(1) As in the previous case, for T > T), the only physical state is given by the 
gaseous phase po, corresponding to zero deviation shown in the figure. 

(2) As the temperature T is lowered crossing Tq, the gaseous state losses its 
stability, leading to two metastable states: one is the liquid phase p + p + , 
and the other is the gaseous phase p +p~ . Again, the gaseous phase p$+p~ 
does not appear, and the asymmetry principle of fluctuations is valid in this 
situation as well. 

(3) The phase transition here is of the second order, as the energy is continuous 
at Tq. In fact, the energy for the transition liquid state is given by 

9 

G(p + P + ) = G(p ) -£-(T- T c f for T < T c , 
4a 3 

for some a > 0. Hence the difference of the heat capacity at T = T c is 

d 2 a 2 

AC = -T c ^ (G($+(T+)) - G(p )) = —T c > 0. 

Namely the heat capacity has a finite jump at T = T c , therefore the tran- 
sition at T = T c is of the second order. 

Fourth, on the curve BC, Oa(T,p) < 0, and the phase transition diagram is 



given by Figure 4.4 



(1) For T > Ti, the system is in the gaseous phase, which is stable. 

(2) For T) < T < Ti, there are two metastable gaseous states given by po and 
Po + P~ ■ As before, although it is metastable, the gaseous state p + p~ 
does not appear. 
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(a) (b) 

Figure 4.3. Continuous transition for the case where a 2 = 0. 

(3) For T < T , the gaseous phase po losses its stability, and the system under- 
goes a dynamic transition to the metastable liquid state po + p + . Mathe- 
matically, the gaseous state po + P~ is also metastable. However, it does 
not appear either, due to the asymmetry principle of fluctuations. 

(4) The dynamic transition in this case is of the third order. In fact, we have 



G(po + P + ) = G(po(T)) - 



6oa 



r (T - T) 3 + o(|T - T| 3 ) for T < T . 



Namely, the free energy is continuously differentiable up to the second order 
at T — To, and the transition is of the third order. It implies that as 
(Tq,pq) £ CB the third-order transition at (T ,po) can n °t be observed by 
physical experiments. 




FIGURE 4.4. Type-Ill (mixed) transition for a 2 < 0. 



In summary, we have obtained a precise characterization of the phase transition 
behavior near the Andrews point, and have derived precisely the nature of the 
Andrews critical point. In particular, we have shown the following: 

(1) The transition is first order before the critical point, second-order at the 
critical point, and third order after the critical point. 

(2) The curve X(T,p) = always defines the gas-liquid co-existence curve in 
both sides of the critical point C. We note that in the classical theory, the 
co-existence curve terminates at the critical point, and with our theory, we 
are able to determine the gas-liquid transition behavior and the co-existence 
curve beyond the critical point. 
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5. Asymmetry Principle of Fluctuations and the Preferred 
Transition Mechanism 

Wc have shown in the last section that in all three cases (both sides of the critical 
point and at the critical point), the metastable state po+p~ does not appear. Hence 
the only possible physical explanation is the asymmetry of the fluctuations. In fact, 
for the ferromagnetic systems we also see this asymmetry of fluctuations 2J . This 
observation leads to the following important principle: 

Physical Principle (Asymmetry of Fluctuations) . The symmetry of fluctuations 
for general thermodynamic systems may not be universally true. In other words, 
in some systems with multi- equilibrium metastable states, the fluctuations near a 
critical point occur only in one basin of attraction of some equilibrium states, which 
are the ones that can be physically observed. 

An alternate explanation of this principle is related to phase transitions in cer- 
tain preferred direction in a given thermodynamic system, which we call preferred 
transition mechanism. We conjecture that this mechanism is universal as well. 
Here we use this mechanism to explain the asymmetry principle of fluctuations in 
the gas-liquid transition, and, in return, to explain the meaning of the preferred 
transition mechanism. 

In the gas-liquid transition as the temperature is lowered, the system prefers 
phase transitions to denser phase. This can be considered as one aspect of the 
preferred transition mechanism. 




Figure 5.1. Preferred transition mechanism for a 2 > with fixed pressure. 
Another important aspect of the mechanism is the preferred transition at a 



critical point T* as shown in Figure 5.1 which is reproduced from Figure |4~2[a) 



This critical point is between T and Ti, and for water under one atmospheric 
pressure, T* is 100 °C. For T* < T < T\, the liquid state state p + p + is called 
superheated liquid, and for T <T < T* , the gaseous state po is called supercooled 
gas. The preferred transition mechanism consists of the following: 

(1) As temperature decreases, the system is forced to undergo a first-order 
transition at T — T* , from the gas state po to the liquid state po + P + ■ 

(2) As the temperature increases, the system is forced to undergo a first-order 
liquid po + P + to gas po transition as the same critical point T = T* . 

(3) The supercooled gas and superheated liquid can rarely occur. Physically, 
this is related to the hysteresis phenomena. 
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